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ABSTRACT 

We present near-IR VLT/ISAAC and mid-IR Spitzer /IRS spectroscopy of the young 
massive cluster in the W31 star-forming region. f/-band spectroscopy provides refined 
classifications for four cluster members O stars with respect to Blum et al. In addition, 
photosphcric features are detected in the massive Young Stellar Object (mYSO) #26. 
Spectroscopy permits estimates of stellar temperatures and masses, from which a clus- 
ter age of ^0.6 Myr and distance of 3.3 kpc are obtained, in excellent agreement with 
Blum et al. IRS spectroscopy reveals mid-infrared fine structure line fluxes of [Neii- 
m] and [Sni-iv] for four O stars and five mYSOs. In common with previous studies, 
stellar temperatures of individual stars are severely underestimated from the observed 
ratios of fine-structure lines, despite the use of contemporary stellar atmosphere and 
photoionization models. We construct empirical temperature calibrations based upon 
the W31 cluster stars of known spectral type, supplemented by two inner Milky Way 
ultracompact (UC) Hn regions whose ionizing star properties are established. Calibra- 
tions involving [Neiil] 15.5//m/[Neii] 12.8/xm, [Siv] 10.5/xm/[Neil] 12.8/xm or [Ann] 
9.0/im/[Ncii] 12.8/xm have application in deducing the spectral types of early- to mid- 
O stars for other inner Milky Way compact and UC Hn regions. Finally, evolutionary 
phases and timescales for the massive stellar content in W31 are discussed, due to the 
presence of numerous young massive stars at different formation phases in a 'coeval' 
cluster. 

Key words: (Galaxy:) open clusters and associations: individual: W31 (G10.2-0.3) 
- (ISM:) HII regions - Stars: early-type; - Stars: fundamental parameters - Infrared: 
ISM 



1 INTRODUCTION 

The formation o f high mass stars remains an unsolved astro- 
phys ical puzzle l|Zinnecker fc Yorkc 2007; C larke fc Bonnefil 
l200ct ). Unlike the situation for low mass stars, for 
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which multiwavelength observations are plentiful, very 
massive stars (A/* >25M©) are born within compact, 
deeply embedded star-forming regions, severely restrict- 
ing observations to either the far-infrared where heated 
dust dominates their appearance, or radio wavelengths 
where gas ionized by the central star(s) can be charac- 
terised. Either way, high mass stars themselves cannot 
be see n directly within these ultracompact (UC) Hn re- 
gions (|Churchwelll l2002h until the column density of dust 
along our line of sight falls below 2-4 magnitudes in 
the _R"-band. To date, on ly a few such cases have been 
identified - G23.96+0.15 (|Hanson. Luhman fc Rieke] 120021 : 
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Figure 1. (a) i^-band image of W31 from lBlum et al] (|200lf) . showing the 4 slit positions for the near-IR ISAAC spectroscopy (dashed 
lines) and the 9 IRS apertures for the first nod position of the mid-IR spectroscopy. North is up, East to the left and the image is 72" 
X 72" or 1.2 X 1.2 pc at our adopted distance of 3.3 kpc; (b) 3.6^m IRAC image of W31 showing both IRS nod positions (first solid; 
second dotted) at the same scale and orientation as for (a). 



Crowther fc Furnesj2008h . G29. 96-0.02 dWatson fc Hanson 
1997| : lHanson et alj|2005h, G45.45+ 0.06 (|Blum fc McGregor 
2008T ), W51d jBarbosa et all 120081 1 - with the stellar con- 



tent of other UC Hn regions r eliant upon indirect far-IR 
or radio continuu m techniques (|Wood fc Churchwelll [l989l ; 
iKurtz et al]|l994 >. 

Fortunately, the advent of efficient mid-infrared imag- 
ing and spec troscopy from space with Infrared Space Obser- 
vator y ( ISO. [Kessler et"al]|l996l ) and Spitzer (|Werner et al] 
120041 ), plus ground-based 8-10m telescopes has opened up a 
further window to study such embedded regions. Specifically, 
a number of fine-structure lines from ionized regions are seen 
in the mid-infrared, notably [Nell-in], [Siii-iv], which al- 
beit indirect, provide information upon the hardness of the 
extreme ultraviolet (EUV) radiation from their constituent 
O stars. Such diagnostics, analogous to the optical forbid- 
den lines of [Oii-iii] and [Sii-iii], may allow the 'inverse 
problem' of establishing properties of the i onizing stars of 
ultra -compact and compact Hn regions (see lOkamoto et al] 
120031 ). Alternatively, indirect approaches based upon near- 
IR hydrogen and helium nebular lines are also e mployed (e.g. 
iLumsden et~aT]|2003l ; iBlum fc McGregorll2009T l. 

To date, tests of photoionization and stellar atmosphere 
models using these mid-IR line diagnostics have been rather 
unsatisfactory. Firstly, this is because Hn regions are usu- 
ally ionized by multiple early-type stars in compact clusters, 
and secondly the strength of mid-IR fine structure lines is 
affected both by the energetic photons from OB stars and 
nebular properties. The most comprehensive study of a UC 
Hn re gion attempted to date has been by iMorisset et al] 
|2002T ) for G29.96-0.02 whose results suggested a tempera- 
ture of T c ff ~ 35 ± 3kK for the ionizing star, in contrast to 
?cft ~ 41 ± 2kK from a non-LTE analysis of near-IR spec- 
troscopy for the star |Hanson. Puis fc Repolu st 2005). 



Photoionization models are widely used to infer the stel- 
lar pr o perties of both embe d ded H 11 regions ((Sellmaier et al] 
19961; iGiveon et al] [20021; |Sim6n-Dfaz fc Stasiriskal 1200a 



Perez-Montero fc VflchezTbOOgl) and entire galaxies from 



mid-I R fine structure lines ( Lutz et alii 19961 ; iRigbv fc Riekd 
|2004| ). If the highly discrepant effective temperatures ob- 
tained for G29.96-0.02 using indirect techniques were re- 
peated for other, single, embedded O stars, then previ- 
ously published results from mid-IR diagnostics may be 
called into question, especially those at high metallicity (e.g. 
iThornlev et alj|200Ch . 

The focus of the present study is the Galactic Giant HII 
region (GHR) W31, specifically th e young star cluster (10.2- 
0.3, h ereafter W31) discussed by Blum. Damineli fc Contil 
|200ll '). This cluster, for which IBlum et all (|200ll ) estab- 
lished a distance of 3.4 kpc and extinction of Ak = 1.7 mag, 
hosts a minimum of four 'naked' O-type stars, as deduced 
from near-IR spectroscopy, plus a number of massive stars 
which are still embedded at near-IR wavelengt hs (hereafter 
mass ive YSOs) plus numerous UC H 11 regions (I Ghosh et al] 
Il989l ). As such, this cluster provides an excellent oppor- 
tunity to study the different early evolutionary phases of 
massive star formation. In addition, it possesses an unusual 
morphology, in that the highest mass stars are located at the 
periphery of the cluster, in contrast to the m ass segregated 
morphology of most ot her young clusters (|de Griis et al] 
120021 ; lAllison et al]|2009h . We shall exploit this unusual ge- 
ometry through spectroscopy of individual ear ly-type stars 
in W 31 with the Infrared Spectrograph(IRS, iHouck et al] 
|2004| ) aboard Spitzer, supplemented with H and K band 
VLT observations w ith the Infrared Spectrog raph and Ar- 
ray CameraQSAAC. lMoorwood fc Cubvlll99sh . 

Our aim is to compare photoionization model predic- 
tions of individual Hn regions with stellar temperatures 
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Figure 2. H band (left panel) and K band (right panel) VLT/ISAAC spectra of the O stars in W31, plus template stars from 
lHanson et al.l l|2005t ). All stars show stellar absorption of Hen and He I in the H band, the ratio of which we use to spectral type the 
stars. The stars also show C IV and N III emission in the i-sT-band indicating early— mid type O stars. The later of the O stars also show 
Br7 emission, indicating a younger evolutionary phase where more nebula material is present around the star. 



obtained directly from near-IR spectroscopic analysis. We 
are thus able to significantly increase the statistics of H n 
regions for which mid-IR diagnostics are available, whose 
ionizing stars are well constrained. Our results are of poten- 
tial application in deducing ionizing stars of UC H II regions 
within the inner Milky Way based upon fine structure diag- 
nostics available from ground -based telescopes , namely [S iv] 
10.5/xm and [Nell] 12.8/nm (|Zhu et al.ll2008l ). The diffrac- 
tion limit of 8-10m telescopes in the N-band is ~0.4 arcsec, 
equivalent to a physical scale of 0.01 pc at a distance of 5 
kpc, well below the .1 pc size of typical UC Hn regions 
ijOkamoto et al.ll2003l ). 

The present paper is structured as follows. § [2] detail 
the VLT and Spitzer observations, while fj3] presents refined 
spectral classification of naked O stars in W31. §[4]presents 
the mid-IR fine structure lines from W31 sources, including 
predictions from photoionization models and empirical cal- 
ibrations. A discussion of mid-IR fine structure lines is pre- 
sented §[5] together with the massive star content of W31. 
Brief conclusions are drawn in § [6] 



2 OBSERVATIONS AND DATA REDUCTION 

Two main observational datasets are used in the present 
study, obtained with the VLT ISAAC near-IR spectrograph 
and Spitzer IRS mid-IR spectrograph. 



2.1 VLT ISAAC spectroscopy 

Long-slit H- and if-band near-infrared spectroscopy of 
sources in W31 were obtained with the ISAAC instrument 



mounted at the Very Large Telescope between 4th April- 
17th June 2006 (Programme 077.C-0550(A), P.I. Crowther). 
Four slit positions were used, as illustrated in Fig.[]Ja). Each 
included two or more targets of interest, namely th e naked 
O stars and massive YSO's from lBlum et alj |200ll ). 



The detector was the 1024x1024 Hawaii Rockwell 
array, while three medium resolution grating settings 
(0.775A/pixel) were obtained, centred at 1.71, 2.09 and 
2.20 /im. These observations were obtained at low airmass 
during variable seeing conditions using a 0.6 arcsec wide 
slit, and reduced using standard iraf packages. We observed 
using ABBA nod-cycles, a standard infrared A number of 
AB pairs were obtained for each grating setting with wave- 
length solutions achieved from comparison XeAr arc images. 
From these, the observations covered 1.671-1.751 /im, 2.029- 
2.155 fim and 2.140-2.265 (im at spectral resolutions of 3.8A, 
6.0A and 6.0A respectively, as measured from arc lines. 



Telluric correction was achieved by spectroscopy of 
early-G dwarfs observed at similar airmass to W31, cor- 
rected for their spectral features using high resolution ob- 
servations of the Sun, adjusted for the radial velocity and 
spectral resolution of the template stars. An extensive dis- 
cussion of telluric correction for medium reso lution near-IR 
spectr oscopy of early- type stars is provided bv lHanson et al.l 
(2005). The 2.09^im setup suffered from low-level variable 
structure which was accentuated upon flat-fielding, arising 
from a 50 Hz pickup inherent to the instrument. Only the 
two other settings were flat-fielded. Consequently, the con- 
tinuum S/N achieved was lower for this grating position 
(~100) than the other settings (~200). 
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Tabl e 1. Classification of early-type stars in W31 from VLT/ISAAC spectroscopy based upon the near-IR calibration of 
ICrowther fc Furnesd feOOSft . 





EW (A) 

Hen (1.692/wn) Hei(1.700/im) 


WMcII 
8 Ho I 


Spectral Type 
This Work Blum et al. 


Further Comments 


#2 
#3 
#1 
#5 


0.59 ± 0.08 
0.40 ± 0.08 
0.58 ± 0.06 
0.68 ± 0.05 


< 0.2 
0.27 ± 0.07 
0.87 ± 0.06 
0.96 ± 0.04 


>0.47±0.06 
0.17±0.20 
-0.17±0.07 
-0.15±0.05 


03-4V 
05±?' 5 V 
O5.5±0.5V 
O5.5±0.5V 


05.5 V 
05.5 V 
05.5V 
05.5V 


N III, C IV in emission 

N III, C IV in emission 

Nm, Civ in emission 

N III, C IV in emission, nebular 

emission from Br7, He I 2.058/xm 



2.2 Spitzer IRS spectroscopy 

W31 was observed with the mid-IR spectrograph IRS in GO 
#3337 (W31CLUST, P.I. Crowther) between 14-16 Septem- 
ber 2005 using all four modules, sampling the short wave- 
length region at low and high resolution (SL and SH), plus 
the long wavelength region at low and high resolution (LL 
and LH). I n addition, W31 was imaged with the IRAC 
instrument (|Fazio et al.l 120041 ') at 3.6, 4.8, 5.8 and 8.0/im, 
using HDR mode, with 12 s exposures in a 12 position, 
Reuleaux dither pattern. Unfortunately, the IRAC 5.8 and 
8.0/im imaging and long wavelength IRS observations were 
heavily saturated, so our analysis largely focuses upon the 
SH, staring mode (30 sec, 6 cycles) observation with IRS, 
obtained at two nod positions. The spectral range covered 
was 9.9-19.6^m, at a resolution of R~600. 

The individual 4.7x11.3 arcsec 2 (2x5 pix) apertures 
from the first nod position are superimposed upon a 
A"— band image of W31 in the FigQJa). FigQJb) shows both 
IRS nod positions together with the IRAC 3.6/im image of 
W31. A background subtraction was applied to these two 
dimensional IRS datasets - prior to extraction - using a ded- 
icated OFF position at 18:09:30.4 -20:20:00 (J2000) from a 
second set of cluster positions which are not presented in this 
paper. This pointing gives a correction for nebular emission 
in the extended Hn region, including the PAH's. However, 
it does not correct for emission from the dense material close 
to each of the targeted stars. 

Observations have been processed and calibrated in the 
standard IRS pipeline to basic calibrated data products 
(BCDs), using a version which is compatible with S16.2. 
Multiple exposures at the same position have been co-added 
at the BCD (2 dimensional) level, and the OFF position 
has been used to correct rogue pixels and subtract the ex- 
tended background. Then the spectra were extracted with 
the offline S16.2 post-BCD pipeline. Spectra from the two 
staring-mode nods were averaged, except for cases for which 
a single nod was used to minimise contamination. Specifi- 
cally, only nod 1 was used for #1 and #26 and only nod 2 
was used for #2 and #30. This is to avoid contamination 
from potentially interloping sources during extraction of the 
spectra (recall Fig.[TJb)). These sources are identified within 
the nods from any of the available imaging observations, in- 
cluding unsaturated regions of the IRAC datasets. 

Since the width of the full width at half maximum of 
the central Airy pattern of the IRS point source function 
is comparable to the size of the SH slit itself, the full slit 
is extracted, so it is not possible to extract (or even easily 
identify) separate sources or the background. 



Each of the echelle spectral orders were trimmed at the 
red edges, corresponding to the lower part of the array where 
photometric response sharply drops. Final spectra have been 
created by merging the orders into single spectra, with ex- 
cellent inter-order agreement when the point-source flux cal- 
ibration is applied, with only minor offsets (well within the 
photometric uncertainties of 15%), indicating that the dom- 
inant source signal from each of the intended targets behave 
more consistently as point sources than extended sources. 
Formally, a la absolute flux calibration of 10% was achieved, 
except for the raw 18.7/im [S ill] line which was saturated 
in a few cases. These were corrected by linearization and 
extrapolation, from which an additional 5% uncertainty re- 
sulted. 

Finally, we should note that these observations have 
been obtained when the SH array was relatively undam- 
aged by the cumulative effects of space weather, thus receive 
good correction for rogue pixels with data at the OFF po- 
sition. Furthermore, these observations are not affected in 
the calibrations or data quality by more recent versions of 
the pipelines (through S18.x) which mainly benefit the IRS 
peak up imaging AOT and data obtained with more aged 
arrays. 



3 PROPERTIES OF EARLY-TYPE STARS IN 
W31 FROM NEAR-IR SPECTROSCOPY 

3.1 Near-IR classification of W31 O stars 

In Figure [2] we present our H- and K-h&nA ISAAC spec- 
troscopy of naked O stars in W31, together with spec- 
tra of templa t e O s t ars from the high re s olutio n atlas of 
lHanson et all (|2005h . ICrowther fc Furnesd (|2008l ) have re- 
cently presented a method of determining O dwarf subtypes 
from the ratio of the observed Hell 1.692/im to He I 1.700/im 
equivalent widths. These were measured from the emission 
line fitting (ELF) suite of routines in the STARLINK spectro- 
scopic analysis programme DIPSO and are presented in Table 
[T] together with inferred spectral types. 

Stars #3, #4 and #5 all show clear He I and Hell 
absorption features, from which 0515' 5 V, O5.5±0.5V and 
O5 .5±0.5V classifications are o btained, according to Fig. 3 
of ICrowther fc Furnessl (|2008T l . in close agreement with 
iBlum et al.l (|200lh . Jf-band spectroscopy reveals C IV, N III 
emission features, plus Br/y and Hen 2.189^im absorption, 
the core of the former filled-in by nebular Br7 emission. 

In the case of #2, negligible He I 1.700/^m is observed, 
from which an 03-4V classification is inferred. If we had 
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Table 2. Interstellar extinctions, Aj$, and distance modulus 
(DM) towards individua l stars in W31 using near-IR photomc- 
tr v of iBlum et al.l d200ll) and the absolute magnitude calibration 
of iMartins fe Pled 12006T) 



ZAMS 10 5 ' 6 10 5 - 8 1 !)" 



W31 m K 


H - K 


(H-K) ( 




A K 


M K 


DM 


# mag 


mag 


mag 


mag 


mag 


mag 




2 10.02 


1.16 


-0.10 


1.26 


2.29 


-4.98 


12.71 


3 10.30 


1.06 


-0.10 


1.16 


2.11 


-4.39 


12.58 


4 10.34 


0.88 


-0.10 


0.98 


1.78 


-4.27 


12.83 


5 10.37 


1.26 


-0.10 


1.36 


2.47 


-4.27 


12.17 


Average 








2.16 




12.57 


(7 








0.29 




0.29 




Wavelength 



Figure 3. //-band and ZC-band spectroscopic fits (dotted lines) 
to VLT/ISAAC (solid lines) observations of hydrogen and helium 
lines in W31 #2 (03-4 V), #3 (05+°- 5 V), #4 (O5.5±0.5 V) and 
#5 (O5.5±0.5V). 



classified #2 on the basis of its K-band spectrum, the pres- 
ence of significant C IV emiss ion would have su ggested a 
spectral type of 04V or later (|Hanson et al.ll2~005l ). 



3.2 Distance and Age of W31 

Armed with our refined s pectral type s , we now use the 



Armed with our refined s pectral type s , we now use the 
near-IR photometry from IBlum et al. I (120011) and obser- 
vational absolute magnitude calibration of IMartins fc Ple3 
|2006t ) to obtain a revised distance to W31. This is pre- 
sented in Table [2] in which observed colours of O stars pro- 



Table 3. Stellar properties of the W31 O stars from non-LTE 
CMFGEN models using a distance of 3.3 kpc. Uncertainties in 
stellar temperatures are ±1.5 kK, with the exception of #2 for 
which an upper limit is 48 kK. current stellar masses are ob- 
tained from co mpar i son w ith 2xZq theoretical isochrones from 
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Schaererl ll200ll ). 
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Figure 4. H crtzsprung-Russell di agram of the 4 naked O stars 
in W31. lLeieune fc Sc hacrcr (200j]) isochrones are shown for ages 
of 10 3 , IP 56 10 5 -8 , 10 6 and 10 6 2 yr, based on the 2xZ e 
iMevnet et al. I l ll994h evolutionary tracks, from which a cluster 
age of ~0.6 Myr is obtained. 



vide a direct measurem ent of interstellar extinc tion from 
A K ~ 1&2E H -K from llndebetouw et all (|2005l ). For #4 
and #5 uncertainties in subtypes have little influence upon 
absolute magnitudes, while subtype uncertainties formally 
introduce absolute magnitude uncertainties of up to 0.3 mag 
for #2 and #3. These are mitigated somewhat by the typi- 
cal spread of ±0.5 mag in absolute magnitudes for individual 
spectral types. 

Overall, we obtain a higher interstell ar extinction of Ak 
= 2.16 ±0.29 mag towards W31 than IBlum et all (|200ll ) 
on the basis of updated intrinsic colours, but this is largely 
cancelled out by the revised absolute magnitude calibration. 
As such, we obtain a very similar overa ll distance (3.3 1° \ 
kpc) to the 3.4 kpc distance obtained bv lBlum et all (|200ll ). 

For an adopted Galactic Centre distance of 8.0 kpc 
(|Reidlll993l ) for the Sun, our preferred distance to W31 sug- 
gests that it lies ~4.8 kpc from the Galactic Centre. The 
Galactic oxygen metallici ty gradient is A log (O/H) = - 
0.044 ± 0.010 dexkpc -1 jEsteban et all 120051 ), from which 
we anticipate that W31 is 40% more metal-rich than Hn 
regions within the Solar circle, i.e. log O/H ±12 ~ 8.81, 
compared to 8.65 for the Orion Nebula at a Galactocentric 
dista nce of ~8.4 kpc (|Esteban et al . 2005; S imon- Diaz et af] 
120061 ) . We are unable to derive a neon or sulphur abundance 
for W31 from the IRS spectroscopy since the nebular elec- 
tron temperature is unknown. 

Having established a slightly refined cluster distance, 
we now estimate stellar temperatures (and in turn lumi- 
nosities) from spectroscopic fits to the near-IR diagnostic s 
with the non-LTE CMFGEN code of lHillier fc Miller] i|l998f ). 
CMFGEN solves the radiative transfer equation in the co- 
moving frame, under the additional constraint of statistical 
equilibrium. Since CMFGEN does not solve the momentum 
equation, a density or velocity structure is required. For the 
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Figure 5. H band (top panel) and K band (lower panel) 
VLT/ISAAC spectra of the massive YSOs in W31. No H band 
spectra are available for W31 #1 and #9. Br7 is seen in emis- 
sion throughout, while underlying photospheric absorption is also 
seen in #26, demonstrating an evolutionary phase between the 
embedded and naked O star phases. 



supersonic part, the velocity is parameterized with a clas- 
sical /3-type law, with an exponent of /3=1 adopted. This 
is connected to a hydrostatic density structure at depth, 
such that the velocity and velocity gradient match at the 
interface. The subsonic velocity structure is set by a corre- 
sponding f ully line-blanketed pla ne-parallel TLUSTY model 
(v.200, see lLanz fc Hubenvll2003l). The atomic model is sim- 
ilar to that adopted in Iffillier et all (|2003l ), including ions 
from H, He, C, N, O, Ne, Si, S, Ar, Ca and Fe, with metal 
abundances inc reased by a fac t or of two relative to those 
of the Sun from I Asplund et al.1 i|2004h for O, Ne and Ar, or 
ICoxl (|2000l ) otherwise. 

We have assumed a depth-independent Doppler profile 
for all lines when solving for the atmospheric structure in the 
co-moving frame, while in the final calculation of the emer- 
gent spectrum in the observer's frame, we have adopted a 
uniform turbulence of 50 kms -1 . Incoherent electron scat- 
tering and Stark broadening for hydrogen and helium lines 
are adopted. Finally, we convolve our synthetic spectrum 
with a rotational broadening profile. In view of the quality 
of our observations, rotational velocities should be reliable 
to ±50 kms" 1 . 

For an adopted (uniform) terminal wind velocity of 2000 



kms^ 1 , surface gravity of log(? = 4 and abundance ratio 
of He/H=0.1 by number, we varied the stellar radius and 
(non-clumped) mass-loss rate until an acceptable match to 
the Hen 1.692/tm, Hel 1.700/im, Br7 and Hell 2.189/imwas 
achieved. For the case of #2 (03-4V), for which negligible 
Hel 1.700/im absorption is detected, we favour a tempera- 
ture of ~45 kK since higher temperatures produce too weak 
Br7 absorption. Spectroscopic fits are presented in Fig. [3] 
while resulting physical and wind properties are shown in 
Table [3] Fits to diagnostic lines are generally satisfactory 
with the exception of #3 for which Hell 2.189/tm suffers 
from low S/N, such that we rely upon the weaker 1.692/im 
line in this instance. Note that prominent nebular emission 
features are seen in the hydrogen Brackett series, with the 
possible exception of #2. 

Our spectroscopic analysis did not yield a precise mea- 
surement of surface gravities (and in turn masses), so hav- 
ing established temperatures and lu minosities, we over- 
lot 2 xZp) theoretical isochrones from iLeieune fc Schaereil 
200 ll ) in Figure [4j from which an age of ~0.6 Myr was 
obtained, together with corresponding surface gravities and 
mass estimates (see Table [3}. Three of the four O stars lie 
along a common isochrone of ~0.5 Myr, although from Fig. [4] 
#5 suggests a greater age of >1 Myr. However, it is possible 
that this source is a close double with equal mass compo- 
nents, whose absolute magnitudes would have Mk = -3.9 
mag for a distance of 3.3 kpc. In this case, individual com- 
ponents would have properties similar to those of #4 in Ta- 
ble [3] reducing the inferred cluster age to ~0.5 Myr. 



3.3 Near-IR spectroscopy of massive YSOs 

In addition to the naked O stars, we have obtained ISAAC 
spectroscopy of four massive YSOs in W31 - high mass stars 
whose photospheric features are veiled by circumstellar dust 
with IR exces ses due to hot dust . We observed #1, #9, #26 
and #30 from lBlum et al] (|200ll ) at both grating settings in 
the K-h&nA spectra, but only #26 and #30 were observed 
in the .ff-band. These datasets are shown in Figure \S\ In 
ge neral, our observat ions confirm the lower quality datasets 
of iBlum et all i|200ll ), with nebular Brackett emission lines, 
plus nebular Hel 1.700, 2.058/mi for #30. 

Uniquely among the massive YSOs, #26 shows clear 
photospheric Br7 absorption in addition to weak nebular 
Br7 emission, plus Hell 2.189/im, 1.692/im and Hel 1.700/im 
absorption. From the observed H-band classification diag- 
nostics in the H-band we estimate an approximate spectral 
type of 06, since log W x (Hell 1.692)/W A (Hel 1.700) ~ 
-0.26. This star appears to be at the point of revealing its 
photosphere, despite the presence of a stron g IR dust excess 
(H-K = 2.72 mag from lBlum et al.ll200lT ). 

In principle, we could estimate the physical properties 
of #26, by combining an estimate of its stellar temperature 
~36 k K) with the ff-band photometry from IBlum et al] 
20011 ). the mean Jf-band extinction for the naked O stars 
(Ah ~ 3.4 mag) and the adopted distance of 3.3 kpc. For this 
star the Jf-band is selected on the basis of significant dust 
emission in the if-band. Unfortunately, the resulting param- 
eters for #26 are unphysical since log L/Lq = 4.2 would be 
obtained, i.e. mass of only 4M©(!) for an assumed surface 
gravity of logg — 4.0. More plausibly, #26 may possess a 
small dust excess, i.e. a large colour excess due to extreme 
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Figure 6. IRS short high (SH) medium resolution spectroscopy of W31 sources, for which the fine structure lines of [Siv], [Nell], [Ncin] 
and [S III] are indicated in #2, and the PAH (UIR) features are marked in #9. These mid-IR datasets have been background subtracted 
using an OFF pointing. 



extinction. Recall that this source lies close to the redden ing 
line in the colour-magnitude diagram (|Blum et al.ll200lh . If 
the (H — K) color were to arise solely from line-of-sight ex- 
tinction, An ~ 7.9 mag, from which log L/Lq — 6.2 would 
be obtained. In reality, the physical properties of #26 will 
lie between these two extreme cases. Regardless, this star 
warrants further study since it uniquely displays both stel- 
lar and circumstellar features amongst early-type stars in 
W31. 



4 INDIRECT PROPERTIES OF EARLY-TYPE 
STARS IN W31 FROM IRS SPECTROSCOPY 

4.1 Mid-IR fine structure line ratios 

We present the Spitzer IRS short high medium resolution 
observations of the naked O stars and massive YSOs in 
Fig [6] Although these datasets are background subtracted, 
no allowance for extended nebular emission close to each 
source is made. Their continua are dominated by warm dust 
emission, while the solid state Polycyclic Aromatic Hydro- 
carbon (PAH) features at 11.3pm, 13.6pm and 16.5pm - 
previously known as Unidentified Infrared (UIR) emission 
- are also seen (the 16.5pm PAH feature is not detected 
in #2). There is also a tentative detection of 15. 2pm ab- 



sorption by icy CO 2 mantles on the silicate dust grains 
(|van Dishoeck|[2~004r i, Shallow profiles with a soft blue wing 
starti ng near 14.9pm a re reminiscent of the W33A massive 
YSO (|Gibb et alj|2004h . although we have no confirmation 
of the associated 4.3pm CO2 absorption which is strong in 
massive protostars. 

We have measured the observed fluxes of the fine struc- 
ture lines in all IRS sources, using the ELF routine in 
DIPSO, which are presented in Table [4] Other, weaker fea- 
tures are observed in the SH module, including Hu a (H I 
7-6) at 12.37pnfl for which log[J(Hu a) //([Nell] 12.8pm)] 
= -1.92 ± 0.04. Where appropriate, we include de-reddened 
line intensities in Table [4] obtained from mean /("-band ex- 
tinctions measured previously together with the mid-IR ex- 



1 At the spectral resolution of IRS/SH the 12.37/^m feature is 
actually a blend of Hi (7-6) and (11-8) for which 1(11-8) /I(7- 
6) = 0.123 according t o Case B recombination theory at T e = 
7,500 K and 10 4 cm" 3 iStorev fc Hummer|[l995l) . 
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Table 4. Mid-IR fine structure line fluxes (intensities in bold font), with units of 10 14 Wm 2 , for the O stars and massive YSOs from 
Spitzer/IRS spectroscopy of W31, together with archival 750/SWS observations of selected UC Hll regions from Peeters et al. (2002) 



Line Fluxes (Intensities) 

Object A K [SIV] [Nell] [Nelll] [SIII] log^^T ^llffiT lo ^T[ 









10.5^tm 


12.8/^m 


15.5/mi 


18.7/mi 
















W31 #1 






0.138±0.002 


9.66±0.10 


1.19±0.01 


11.4±0.1 


















2 


1 


0.473±0.006 


12.7±0.1 


1.46±0.02 


13.7±0.1 


-0.94±0.01 


-0 


.46±0 


.01 


-1 


.46±0 


.01 


W31 #2 






0.452±0.004 


7.43±0.07 


3.96±0.03 


12.7±0.1 


















2 


1 


1.55±0.01 


9.78±0.1 


4.86±0.04 


15.4±0.1 


-0.30±0.01 


-1 


00±0 


01 


-0 


80±0 


01 


W31 #3 






0.193±0.002 


5.44±0.05 


1.61±0.02 


8.7±0.1 


















2 


1 


0.66±0.01 


7.2±0.1 


2.0±0.02 


10.5±0.1 


-0.56±0.01 


-1 


20±0 


01 


-1 


03±0 


01 


W31 #4 






0.365±0.003 


7.34±0.07 


3.14±0.04 


11.6±0.1 


















2 


1 


1.25±0.01 


9.66±0.1 


3.85±0.05 


14.0±0.1 


-0.40±0.01 


-1 


05±0 


01 


-0 


89±0 


01 


W31 #5 






O.OlOiO.0001 


5.86±0.7 


1.06±0.01 


6.80±0.1 


















2 


1 


0.341±0.003 


7.71±0.08 


1.30±0.01 


8.2±0.1 


-0.77±0.01 


-1 


38±0 


01 


-1 


35±0 


01 


W31 #9 






0.062±0.001 


7.06±0.11 


0.80±0.01 


8.0±0.1 


















2 


1 


0.211±0.003 


9.3±0.1 


0.99±0.01 


9.7±0.1 


-0.97±0.01 


-1 


66±0 


01 


-1 


64±0 


01 


W31 #15 






0.301±0.002 


8.04±0.08 


3.02±0.03 


12.8±0.1 


















2 


1 


1.03±0.01 


10.6±0.1 


3.71±0.04 


15.4±0.1 


-0.46±0.01 


-1 


18±0 


01 


-1 


01±0 


01 


W31 #26 






0.396±0.004 


13.5±0.1 


3.92±0.03 


15.2±0.2 


















2 


1 


1.36±0.01 


17.8±0.2 


4.81±0.04 


18.3±0.2 


-0.57±0.01 


-1 


13±0 


01 


-1 


12±0 


01 


W31 #30 






0.286±0.002 


12.4±0.2 


3.15±0.02 


14.5±0.1 


















2 


1 


0.98±0.01 


16.3±0.2 


3.87±0.03 


17.5±0.2 


-0.63±0.01 


-1 


25±0 


01 


-1 


22±0 


01 


G23. 96+0.15 






0.33±0.08 


34±1 


2.4±0.2 


25.6±0.4 


















2 





1.3±0.2 


46±1 


3.0±0.3 


31.7±0.5 


-1.19±0.06 


-1 


37±0 


05 


-1 


54±0 


05 


G29.96-0.02 






4.2±0.4 


99±7 


27 ±1 


46±1 


















1 


6 


10.7±1 


122±9 


32±1 


53±1 


-0.58±0.05 


-0 


69±0 


05 


-1 


06±0 


06 



-0.5 



5- -1.0 

CO 



-1.5 
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#26 



G23.96+0.15 ^ 30 #3 
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#4 
#15 



#9 
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Figure 7. Fine structure lines of the O-type stars (filled circles) 
and massive YSOs (open circles) in W31, with several compact 
and ultracompact Hll regions discussed in this work included 
for comparison (open tr i angles) . The dashed line is taken from 
iMartm-Hernandez et al. I ll2002bh . and represents the best fit to 
observational datascts spanning a wide range of environments (see 
text). 



tinction law of I Morris et al.l (|2000t ). namely 
0.637 A K 
0.142 A K 



Of the mid-IR diagnostics, the [S iv] line at 10.51 /im is sig- 
nificantly affected by the 9.7/im silicate feature. The 11.3/im 
PAH emission feature spans a range of equivalent widths, 
typically W\ = 600±300 A, with the exception of two mas- 
sive YSOs, #1 (W x w 2500A) and #9 (Wx ~ 2100A). Weak 
UIR emission a 12.7/xm can be seen in these two sources (see 
Fig. El). 

In addition to the W31 sources, we also present mid- 
IR fine structure line fluxes for selected UC H n regions in 
TableU Of these, G23.96+0.15 and G29.96-0.02 lie within 
the inner Milky Way and possess b oth mid-IR spectros copy 
from Infrared Space Observatory ([Peeters et al] l2002h and 
near-IR spectroscopy of the io nizing star ( Hanson et al.l 
120051 ; ICrowther fc Furnessl 120081). Interstellar e x tincti on for 
these were taken from ICrowther fc Furnessl (|2008h and 
IMartm-Hernandez et al.l |2002al 'l. respectively. 

In Fig. [7] we compare the intensity ratios of 
[Ne ill] /[Nell] and [Slv]/[Sni] for the W31 sources plus se- 
lected compact and UC Hll regions. The former lie close 
to the best fit of observed ratio s in Hll regions in the 



Milky Way and Magellanic Clouds (|Martfn-Hernandez et al 
2002H), bas ed on observat i ons p resented by IPeeters et al 
|2002l ) and IVermeii et ail (|2002r >. The compact and UC 
Hll regions are appa r ent ou tliers, but Fig. 1 from 
IMartm-Hernandez et al.l (|2002b| ) reveals a typical scatter of 
0.2 dex and they did not correct line ratios for interstellar 
extinction. Physically, these are probably amongst the high- 
est density (ionization paramet er) Hll regions of the full 
IMartm-Hernandez et all |2002bl ) sample. 

The W31 sources with the highest ionization are the 
naked O stars #2 (03-4 V) and #4 (O5.5±0.5 V) with mas- 
sive YSOs #1 and #9 possessing the lowest ionization neb- 



r\ 1 0/7 A 
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ulae, with a mixture of naked O stars and massive YSOs at 
intermediate ionization. From these observations, #15, #26 
and #30 could be interpreted as possessing similar stellar 
properties to #3 (Obtl' 5 V) and #5 (O5.5±0.5 V) with later 
subtypes for #1 and #9, although we defer a formal esti- 
mate of their temperatures until § 14.21 Despite this overlap, 
it may be significant that the naked O stars tend towards 
higher ionization nebulae, i.e. higher temperatures (larger 
stellar masses). 

4.2 Photoionization Modelling 

The primary aim of our study is to compare the direct 
near-IR stellar signatures of O stars in W31 (and UC H II 
regions) with the indirect mid-IR nebular lines through 
predictions from photoionization models. We use version 
08.00 of the p h otoion ization code CLOUDY, last described by 
iFerland et al.l (1 19981 '). This solves the equations of thermal 
and statistical equilibrium for a model nebula, represented 
by a sphere of gas with uniform density n and filling factor 
e with a small central cavity which is ionized and heated 
solely by the UV radiation of a single central star. 

Nebular fluxes are predicted, given input abundances, 
ionizing flux distributions and physical parameters, most im- 
portant of which is the ionization parameter 

U = Q /{^Rlnc) 

and Qo is the number of ionizing photons below the H- 
Lyman edge at 912A. Here, Rs is the radius of the Stromgren 
sphere. Alternatively, 

1/3 



u = l 

c 



na B Qg 
36tt 



where ae is the Case B recombination coefficient. For a 
given energy distribution of the ionizing radiation field, any 
combination of parameters which keeps Qu_ne 2 constant will 
result in an identical ioniz ation structure of the gas (see 
IStasiriska fc Leitherei 19961 ). 

Mor isset et al. I~ i|2004 ) have compared a number of stel- 
lar atmosphere codes to mid-IR observations of Milky Way 
Hll regions, conc luding that the non-LT E codes CMFGEN 



2001 ) provide the best 



Simon-Diaz fc S tasiriska 



and WM-Basic (|Pauldrach et al.l 
matc h to observations (see also 

120081 ). We therefore utilise CMFGEN to provide the ion- 
izing flux distributions, as discussed above. 

Ideally, one would employ compact and ultra-compact 
H II regions that were both spherical and ionized by a single 
dominant source for such a study. However, such cases are 
incredibly rare, due to the lack of accurate subtypes for the 
ionizing stars and scarcity of space bourne mid-IR observa- 
tions. G29.96-0.02 satisfies the ideal criteria relative l y wel l, 
but G23.96+0.15 is irregular l|Wood fc Churchwelll 1989h . 



and there is ionized gas throughout the W31 cluster. For- 
tunately, the O stars within W31 are relatively uniform in 
their ionizing output. In contrast to the ONC, where 8 1 Ori 
C dominates the extreme UV radiation field, we assume that 
the ionized gas within IRS apertures centred upon individ- 
ual O stars are dominated by these stars. In reality, the 
diffuse radiation field likely reflects a combination of the 
ionizing photons from multiple cluster members, a consid- 
erations which should be bourne in mind in the following 
analysis. 



4-2.1 Nebular densities 

For the stars whose near-IR spectra enables a spectral type 
to be determined, we can obtain empirical ionization pa- 
rameters if the distance and electron density is known. Ta- 
ble [5] provides ionizing parameters for the two UC H II re- 
gions whose ionizing stars have been determined from near- 
IR spectroscopy, G29.96-0.02 and G23.96+0.15, implying 
log U ~ —1.5 in both cases. 

Specifically, for G29. 96-0.02, lHanson et all (|2005l ) es- 
timated log Qo = 49.6 for a kinematic distance of 7.4 
kpc, which together with a radio-derived Strom gren ra- 
dius of 3.5" (0.13 pc IWood fc Churchwelll Il989l ) requires 
a high electron density of 20,000 cm~ 3 . This is typical 
of UC Hll regions, but is significantly higher than the 
density of 817 cm~ 3 ob tained from ISO observa t ions us - 
ing [Om] 88/52^m by iMartfn-Hernandez et ail (|2002al ). 
For G23.96+0 .15, we adopt the ki n emati c distance of 4.7 
kpc for which ICrowther fc Furnessl (|2008h obtained log Qo 
= 49.4 a nd we adopt a Stromgren radius of 2" (0.025 
pc) from IWood fc Churchwelll (| 19891 ). although this UC 
Hll region has an irregular radio morphology. Again, the 
corresponding electron density of 70,000 cm -3 is much 
higher than the average density of 1543 cm -3 obtained by 
IMartfn-Hernandez et alj (|2002al ) . arising in part due to the 
large aperture of ISO/LWS. 

It was our intention that the IRS/LH spectroscopy of 
W31 would have enabled us to deduce electron densities. Un- 
fortunately, due to severe saturation because of the bright 
dust continuum this was not possible (§ 12. 2p . We have there- 
fore comb ined the radio-derived io nizing flux of logQo = 
50.4 from IConti fc Crowtherl <|2004 ) with the 30" (0.5 pc) 
radius of the Hll region from diffuse Br7 emission (recall 
Fig. [1} to estimate an electron density of 6,500 cm -3 . The 
ionization parameter implied in all cases is close to log U = 
-1.5, albeit slightly lower than for the UC Hll regions, with 
similar valu es anticipated for the massive YSOs in W31. For 
comparison, iBaldwin et~aH i|l99ll ) also derived log U = —1.5 
for the Orion Nebula Cluster (ONC). 



4-2.2 Predicted mid-IR fine structure line ratios 

Table [5] includes predicted mid-IR fine structure line ratios 
for each case using CLOUDY. These wer e obtained with stella r 
models for individual stars from either lHanson et al. I (|2005h . 
ICrowther fc Furnessl (|200ST ) or the present study and were 
calculated with abundances scaled to 1.5 times the Solar 
value, with dust neglected. From a comparison between the 
observed line ratios in Table [4] and predicted values in Ta- 
ble[5]it is apparent that agreement is very poor. For example, 
the [Nem] 15.5/um/[Ne n] 12.8^m and [Siv] 10.5/um/[S in] 
18.7/im ratios differs by 1.5±0.5 dex in all cases. This sug- 
gests either that the stellar or photoionization models are at 
fault, or the problem arises from local background variations 
close to individual W31 sources. The former would appear 
to be more plausible, since the latter explanation would not 
reconcile differences between ISO observations of UC Hll 
regions and predictions. 

Regardless of the origin for the discrepancy, we have 
calculated a grid of photoionization models for three ioniza- 
tion parameters, log U = -1, -2 and -3 using stellar atmo- 
spheric models appropriate for solar metallicity O dwarfs, 
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Table 5. Predicted mid-IR fine structure line ratios for inner Milky Way Hn regions. Elemental abundances are scaled to 1.5 Zq 
throughout. An electron density of 6,500cm -3 is used for W31 on the basis of the radio-derived ionizing flux (lO 50 4 phs" 1 ) and radius 
of ~30" (0.5 pc at a distance of 3.3 kpc) for the giant Hll region. 



Star Sp. Type T cff d logQ n e R s logU log ^"j 1 logiffmj log /' [ Neil] lo S i [Nciij lo Sj\ 

kK kpc cm -3 pc 



W31 #2 


03-4V 


45 


3.3 


49.6 


6,500 


(0.24) 


-1.6 


1.01 


0.42 


0.52 


1.22 


-0.01 


W31 #3 


OS±?- 5 V 


13 


3.3 


49.3 


6,500 


(0.19) 


-1.7 


0.65 


0.23 


0.28 


0.81 


-0.28 


W31 #4 


O5.5±0.5V 


41 


3.3 


49.0 


6,500 


(0.15) 


-1.8 


0.24 


0.01 


0.05 


0.39 


-0.52 


W31 #5 


O5.5±0.5 V 


11 


3.3 


49.3 


6,500 


(0.19) 


-1.7 


0.31 


0.09 


0.08 


0.49 


-0.48 


G29.96-0.02 


04-5V 


11 


7.4 


49.6 


20,000 


0.13 


-1.5 


0.95 


0.60 


0.49 


1.06 


0.00 


G23.96+0.15 


Q7.5V 


38 


4.7 


49.2 


70,000 


0.05 


-1.4 


0.02 


0.46 


-0.06 


0.02 


-0.52 




Figure 8. Log n(S-Ne) versus log [NeIII]/[NeII] for solar metallicity models from [Simon-Dtaz fc St asiriska (2008) (left panel, log U\ n = 
-1, -2, -3) and this work (right panel, log U = -1, -2, -3), together with empirical data for naked O stars (filled circles, identified in left 
panel) and massive YSOs (open circles, identified in right panel) in W31, plus UC Hll regions (triangles, identified in left panel). See 
text for a discussion of the key differences between these photoionization models. 



as listed in Table [6] These were obtained from the O 
star calibration from Con t i et a l. (2008) together with the 
IVink. de Koter fc Lamersi (|200ll ) mass-loss prescription. A 
uniform rotational broadening of 200 kms -1 was applied to 
synthetic spectra. Dust was not included, but has little ef- 
fect upon the predicted line intensities with Orion-like dust 
grain compositions for dust-to-gas mass ratios of 0.1-2%. In 
the unlikely case of Milky Way ISM dust grains with a dust- 
to-gas mass ratio of 2%, line intensity ratios may be affected 
by up to factor of two. 

Observed line intensities are intrinsically linked to 
both ionization parameter, U, and stellar temperature, 
T c ft. One means of isolating the form of the radiation 
field is through the ionization softness parameter r\ = 



([On]/[Om 
optically by 

ness of mid-IR nebular diagnostics can be characterised by 



)/([Sii1/[Siii1) defined for H II regions detected 
Vflchez fc Pagell 1 19881 ). Analogously, the hard- 



r?(S - Ne) = 



([S IV]10.5/j m/[S III]18.7m m) 
([NeIII]15.5/im/[NeII]12.8Atm) 



[Ne ill] /[Neil]. In Fig. [8] we compare our CLOUDY predic- 
tions for Solar composition nebulae for log U — -1, -2, -3 to 
the observed positions of W31 star s and UC H H regions, to- 
gether with recent predictions from lSimon-Dfaz fc Stasiriskal 
l|2008h which are also based on cloudy and CMFGEN mod- 
els. The W31 sources lie very close together in Fig. [8] with 
the marginal exception of the massive YSO #9. This source 
is also an outlier in Fig. [7] possibly as a result of a reduced 
density (ionization parameter). 

G29. 96-0.02 lies close t o the (T eB , log U in ) = (40kK, -2 ) 
photoionization model from lSimon-Dfaz fc Stasiriskal (|200ST ) . 
in reasonable agreement with the empirical properties of 
(41kK, -1.5) from Table [S] In contrast, our results predict 
rather poor agreement for G29.96-0.02 since it sits close to 
the (37kK, -2) p hotoionization model. Th e reaso n for this 
difference is that ISimon-Dfaz fc Stasiriskal (|2008h adopt a 
non-standard model dependent definition of ionization pa- 
rameter, namely 

U in = Q /(4irR? n nc) 



as discussed by Morisset (2004). Variations in U and 
T c ff are broadly independent in a plot of r?(S— Ne) versus 



which is based upon the inner radius of the model cloud, 
Rin for the photoionization calculation rather than the 
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Figure 9. log rj(S-Ne) versus log I[NcIII]/I[NcII] ratio for 2x 

Solar mctallicity grid for logU = —1,-2,-3. Symbols have the 
same meaning as in Fig. 8. 



Stromgren radius, Rs- Consequently, the poor agreement 
for G29.96-0.02 obtained here simp ly confirms the previ- 
ous (poor) consistency discu s sed by iMorisset et al.l (|2002h 
and iMartm-Hernandez et al.l (|2002bl ). The comparison for 
G23.96+0.15 fares little better, with (T eff , log f7) = (36kK, 
-2) predicted, in contrast to empirical estimates of (38kK, 
-1.4) from Table [5] The situation for the W31 O stars is 
fairly similar, with ionization parameters of log U <-3 pre- 
dicted in the right panel of Fig. [5] yet -1.7 is more typical 
of compact and UC H II regions (Table [5} . 

In § 13.21 we indicated that W31 is anticipated to be 
metal-rich with respect to the Solar neighbourhood. Similar 
arguments apply for G23.96+0.15 and especially G29.96- 
0.02. For a kinematic distance of ~7 kpc to G29.96-0.02, 
this UC H II region would lie at Rgc = 4 kpc from the Galac- 
tic Centrfl versus 3.7 kpc for G23.96+0.15 for a kinematic 
distance of 4.7 kpc. 

Therefore, we have also calculated a set of super-solar 
(2 Zq) cloudy models, based upon O star models that 
are identical to the solar composition grid except that both 
metal abundances and mass-loss rates are increased (see Ta- 
ble H5] for the latter). Solely early O stars provide signifi- 
cant numbers of extreme UV photons capable of producing 
[Nem] and [Siv] nebular emission while [Nell] and [S III] 
emission is expected to cease beyond ~09.5V. 

Predictions from 2Zq photoionization models are pre- 
sented in Fig. [9] In general, from a comparison with the Solar 
models, these yield slightly higher stellar temperatures and 
ionization parameters. G29.96-0.02 lies close to (T e g,logU) 
— (38kK, -2), representing a slight improvement with re- 
spect to the empirical temperature. Similar comments apply 
for G23.96+0.15, although predicted temperatures remain 
too low, and ionization parameters are offset by -1.0 to -1.5 
dex for G23.96+0.15 and the O stars in W31. 

2 ISewilo et al.l (120041 ) obtained a distance of 8.9 kpc or Rgc = 
4.5 kpc to G29.96-0.02. 



Table 6. Stellar parameters for solar (twice solar) O dwarf model 
atmospheres used in photoionization calculations, based upon cal- 
ibrations presented in Conti et al. (2008) with mass-loss rates 
obtained from the Vink et al. (2001) prescriptions. 
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Figure 10. Comparison between the predicted [Nelll]/[Nen] in- 
tensity ratio and effective temperature for logl/=-1.5 models at 
Solar and 2 X Solar composition. Empirical data are naked O stars 
in W31 (filled circles) and inner Milky Way UC H II regions (open 
triangles) together with a first order fit to the empirical datascts 
(dotted line, see text). 



4.3 [Ne ill] /[Nell] calibrations 

Let us now focus upon the predicted [Ne lll]/[Ne n] ratio 
versus stellar temperature, which is presented in Fig 1101 for 
the Solar and 2xZq metallicity grids at logf/ = —1.5. We 
have selected this value for the ionization parameter since 
this is characteristic of compact and especially UC Hll re- 
gions in our sample. Empirical results are also shown, with 
stellar temperatures drawn from Table [S] quantifying earlier 
discrepancies with respect to the predictions. An explana- 
tion for the disagreement is not readily apparent, although 
it may involve incomplete line blanketing in current non- 
plane parallel model atmospheres at high energies. Never- 
theless, observations allow us to estimate the temperatures 
of the ionizing O stars in other dense H II regions (for which 
log U ~ —1.5) within the inner Milky Way, to within ±2 kK 
(i.e. one spectral subtype) in principle. A first order fit to 
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Table 7. Estimates of spectral types of the ionizing stars responsible for compact and ultra-compact H II regions in the inner Milky Way 
(-Rgc < 7 kpc) fr om Peeters et al, (2002) plus ma ssive YSOs in W31 from the present study. Interstellar extinctions are taken from 
the present study, [Martm- Hcrnandcz et al. ( 2002a]) , or Ak = 1.6 mag otherwise (marked with *). Kinematic distances adopt a Solar 
galactocentric distance of 8.0 kpc l|Rcid 1993). Temperatures and spectral types are obtained for both the [Ne III] /[Nell] and [Slv]/[Nell] 
diagnostic ratios. 



Source IRAS Alias VLsr d R GC Ref A K log T off Sp Type log T cff Sp Type 

kms" 1 kpc kpc mag ^//^ kK [ n'bVi] kK 
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the empirical datasets using IDL's polyfit routine reveals 

{T ctt /kK) 2 z Q = +45.47+ 1.37 

+6.27+ 1.96(log/[NeIII]//[NeII]). 

O subtypes may then be estimated from the recent T e g 
calibration of lMartins. Schaerer fc Hillierl (|2005l ). Of course, 
this approach is only practical for compact and UC Hn re- 
gions ionized by early- and mid- O stars, in view of the 
weakness of [Nem] 15.5/^m for late O-types. 

To illustrate its potential diagnostic role, in Ta- 
ble [7] we provide stellar temperatures derived from ex- 
tinction corrected [Ne ill] /[Nell] ratios for compact and 
ultra-compact H H regions from ISO/SWS observations 
jPeeters et al] 120021 ). plus the massive YSO's in W31 
from the present study. Ex ti nctions are either taken 
from iMartm-Hernandez et al] (|2002al ). or A K = 1.6 
mag is adopted otherwise. We omit the peculiar com- 
pact Hll region G93. 53+1.47 (Ml-78) from this sample 
|Martm-Hernandez et al] [20081 ). For the massive YSO's in 
W31, Table [7] presents spectral types inferred from the mean 
interstellar extinction obtained from the naked W31 O stars. 
In reality, higher extinctions may be expected for these cases 
(recall § 13. 3[) . Fortunately, the use of a higher extinction does 
not affect the resulting stellar temperatures/subtypes. For 
#26, a subtype of 04-5V is obtained, versus an approximate 
subtype of 06V from its H-band ISAAC spectrum (§ 13.31) . 

Since our calibration has been established using solely 
compact Hll regions within the inner Milky Way - with 
super-solar elemental abundances - one should consider sep- 
arately the case of Hll regions close to, or exterior to, the 
Solar circle. Of course, metal content is very relevant to gas 
cooling and metal-line blanketing of the extreme UV energy 
distributions of O stars. Unfortunately, among the compact 
Hn regions close to the Solar circle, very few possess both 
mid-IR spectroscopy and well determined properties for the 
dominant ionizing star. The ONC - for which 6 1 Ori C (HD 



37022, 06-7 Vp) is the dominant source of ionizing radiation 
- is one s uch case in that it pos sesses both mid-infrared spec- 
troscopv lSimpson et al] (Il998h and a contempor ary analysis 
of its stellar content by ISirnon^ Diaz et all (|2006l ). 

Other examples are generally complicated by the pres- 
ence of multiple ionizing stars of uncertain spectral type 
(e.g. NGC 7538) or mid-IR spectroscopy is lacking. The 
only other case for which mid-IR spectroscopy is available 
together with a well determined spectral type of the ion- 
izing star is Gilt ) . 10+0 .05 (IC 1470, Sh 2-156) for which 
iHunter fc Massevl (|l990l ) derived 06.5 V for its ionizing star. 
Therefore, we refrain from attempting a solar metallicity 
calibration at this time. However, we are able to com- 
pare the ONC with the CLOUDY predictions, recalling that 
iBaldwin et al] (|l99ll ) derived a value of log U = -1.48 for 
the ONC. In contrast to the metal-rich H,ll regions, the ob- 
served neon ratio for the ONC is in better agreement with 
the sol ar-metallicity prediction s for T c g = 39+1 kK for 8 1 
Ori C (|Sim6n-Dfaz et al]|2006l ). 

4.4 [Slv]/[Nen] and [Ar ill] /[Nell] calibrations 

Of course, neither r;(S-Ne) nor [Ne ill] /[Nell] are avail- 
able from ground-based observations. In such cases, only 
[Arm] 8.9pm, [Siv] 10.5/xm and [Neil] 12.8/xm are accessi- 
ble. In Figure [TT1 we compare the [Slv]/[Nen] ratio versus 
stellar temperature for sources in W31, G29.96-0.02 and 
G23.96+0.15, plus CLOUDY predictions for 2 x Solar stellar 
models. 

The comparison between the predicted [Slv]/[Nen] 
ratio and empirical results is much poorer than for the 
[Ne ill] /[Nell] ratio. As such, one might conclude that no 
straightforward means of determining temperatures exists 
from ground-based data alone. However, it has been es- 
tablished that there is a reasonably tight correlation be- 
tween the [Siv]/ [Neil] and [Nem]/[Neii] ratios, as shown 
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Figure 11. Effective temperature versus log [Slv]/[Nell] for our 
log U = —1.5 solar and 2x Solar models together with inner Milky 
Way sources and resulting empirical calibration (dottcd-linc). 

by iGroves. Nefs fc Brandll (|2008h . For the sample of Milky 
Way compact H n regions observed by IPeeters et alj (|2002n 
we obtain a linear, least-squares fit 

log /[Ne III] II [Nell] = 0.86 ± 0.06(log J[S IV] //[Ne II]) 
+0.21 ±0.06 

after correcting line fluxes for extinction (as discussed 
above) . This is p resented in Fig. 1121 For comparison, 
I Groves et al. | (|2008f ) obtained coefficients 0.74 and 0.42 for 
65 Galactic H II regions from lGiveon et all (|2002l ), neglecting 
interstellar extinction. 

We use this conversion to apply our empirical 
[Ne ill] /[Neil] temperature calibration to [S lv]/[Ne n]. This 
is also presented in Fig. 1111 from which similar consistency 
with empirical results to that achieved for [Ne m]/[Ne n] is 
obtained, revealing 

(T ctt /kK) 2 z e = +46.79 ± 1.51 

+5.39 + 1.69(log/[SIV]//[NeII]) 

for the super-solar case. Again, we include the predicted sub- 
types of Galactic compact Hn regions from our empirical 
[Slv]/[Nen] calibrations in Table [7] Since [Siv] is relatively 
sensitive to extinction, if we were to use extinctions inferred 
from their H-K colours (e.g. A K = 4.9 mag for W31 #26) 
rather than the average of the naked O stars, we would ob- 
tain temperatures up to 3kK higher. It is likely that their 
actual extinctions lie between these two extremes due to the 
contribution of hot circumstellar dust. 

In principle, either the [Ne lll]/[Ne n] or [S lv]/[Ne n] cal- 
ibration may be applied to Hn regions ionized by early 
or mid-type O stars, since late-type O stars possess too 
low photon fluxes at the energies required to produce ei- 
ther Ne 2+ (> 41eV) or S 3+ (> 35eV). Therefore, we have 
also considered the potential role of [Arm] 9.0/im from the 
two UC Hn regi ons (this line was no t covered by IRS/SH 
datasets). For the lPeeters et all l)2002h sample of H II regions 
as discussed above, we obtained [Arm] line intensities using 




-2 -1 



log I[SIV]/I[NeII] 




-1.0 -0.5 



log I[ArIII]/I[NeII] 

Figure 12. (Upper panel) Correlation between I[S lv]/7[Nc II] 
and /[Ncm]//[Nell] for Hn regions located in the inner (filled 
squares) and outer (open squares) Milky Way using ISO/SWS 
fluxes from Peeters et al. (2002), dereddened as described in 
the text; (lower panel) As above except for correlation between 
/[Arm]//[Ncn] and /[Nem]//[Nell]. 
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Figure 13. Effective temperature versus log [Ar iii]/[Nc iv] for 
our logf/ = — 1.5 solar and 2x Solar models together with inner 
Milky Way sources and resulting empirical calibration (dotted- 
line). 
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an extinction correction of 
Ag.ofim = 0.719 A K 

and obtained a similar quality fit to that of I[S iv]/I[Nen] 
above (see Fig. 1 12 j) . namely 



log /[Ne III] /I [Nell] 



2.06±0.21(log7[ArIII]/I[NeII]) 
+0.55 ±0.11 



allowing estimates of temperatures to be made from the 
[Arm] /[Nell] ratio. Unfortunately, this is relatively insensi- 
tive to temperature and so should only be considered if nei- 
ther of the primary calibrations are available. Nevertheless, 
the comparison between empirical results for G23.96+0.15 
and G29.96-0.02 and our metal-rich calibration, presented 
in Fig. 1131 is reasonable. For completeness, the super-solar 
calibration may be expressed as: 



{T eS /kK) 2 



+48.92 ± 1.79 

+12.92 ± 3.47(log/[ArIII]/7[NeII]) 



It should be re-emphasised that these calibrations are solely 
intended for compact and ultra-compact H n regions for the 
inner Milky Way, and are liable to revision once additional 
empirical results become available. Nevertheless, it is ap- 
parent that the use of any mid-IR nebular diagnostics to 
establish properties of metal-rich O stars will lead to stellar 
temperatures that are too low, by as much as 5kK or 10%. 



5 DISCUSSION 

5.1 Massive Star Content of W31 

Let us now return to W31 itself, and specifically its massive 
star content. From Fig. [7] the highest ionization stars are 
naked O stars including the highest mass star in W31 #2, 
while the lowest ionization stars (#1, #9) are massive YSOs 
with presumably the lowest masses. 

Ideally, one would like to comment upon differences in 
timescales over which circumstellar material is destroyed by 
young O stars of different masses, but the range of masses 
exhibited by our sample is likely to be relatively modest 
(with the exception of #2) and the unusual geometry of W31 
suggests that massive star formation across this cluster was 
not necessarily coeval. The only naked O star from Fig. [3] 
suggesting an age spread is #5. However, as discussed in 
§ 13.21 it is plausible that #5 is a close binary with similar 
mass components from which a uniform cluster age would 
be obtained. 

For the case of a genuinely coeval cluster, it is likely that 
stars with the highest luminosities - and hardest extreme 
UV photon fluxes - destroy circumstellar dust more rapidly 
than lower luminosity, softer extreme UV stars, which are 
conceivably the massive YSOs. From Table [7] stellar sub- 
types of 04-6V are obtained from our [Nelll]/[Nen] cali- 
bration (Table [7J, i.e. comparable to the naked O stars. 

W31 #26, in particular, appears to be in the process 
of actively clearing its circumstellar dust, revealing photo- 
spheric absorption lines, consistent with a subtype of 06V, 
slightly later than the indirect calibration. #15 and #30, on 
the basis of similar neon ratios to #26 its mid-IR neon lines 
likely possess similar subtypes, with ~06-6.5 subtypes for 



#1 and #9. The lowest (current) mass estimated amongst 
the naked O stars is ~ 36M for the 05.5V star #4 (Ta- 
ble^. Consequently, it is likely that the massive YSO's pos- 
sess stellar masses of order 30-35 Mq. Recall also that W3 1 
also hosts a number of UC H II regions (|Ghosh et alj|l989h . 
Such sources are still too deeply embedded to be detected at 
near-IR wavelengths, most likely as a result of still lower lu- 
minosities/masses (late O subtypes?). If the UC Hn regions 
are genuinely coeval with the naked O stars and massive 
YSO's in W31, a timescale in excess of 0.5 Myr is required 
for late O stars to sufficiently clear their environment to be 
detected at near-IR wavelengths. 



5.2 Analysis of embedded stellar populations 
from mid-IR fine structure lines. 

Our study has focused upon compact and ultra-compact 
Hn regions within the inner Milky Way. We find that all 
mid-IR nebular diagnostics will lead to stellar temperatures 
that are too low, by as much as 5kK or 10%. However, 
mid-IR fine structure lines also have application in exter- 



ral g 


alaxies (IRubin et al. 


2007), starbursts (|Thornlev et al.1 


l200d: IVerma et al.l 20031; 


Brandl et al.l 20061) and ultralumi- 



Within this diverse sample, the mid-IR has the great- 
est diagnostic role for the highly obscured cases, either 
young, embedded massive clusters within metal- poor star- 
bursts such as NG C 5253 (|Crowther et alJll999T). He 2-10 



JVacca et al.1 12002T ). or metal-rich ULIRG's tjCenzel et all 
Il998r i. 

Of course, the aim of studies of starbursts is generally to 
obtain ages and/or stellar content (mass, Initial Mass Func- 
tion) yet similar approaches to the present study are usu- 
ally followed, involving photoionization models and ionizing 
flux distributions from stellar atmosphere models coupled to 
evol utionary predictions through population synthesis codes 
(e.g. iPindao et afll2002ft . 

Problems high li ghted here and elsewhe re 
|Morisset et al.l |2002l : iMartm-Hernandez et all l2002al lbT) 
for Hn regions in the inner Milky Way probably have 
their origins in both (i) incomplete opacities from all 
relevant ions in stellar atmosphere models allowing for 
non-LTE effects and stellar winds; (ii) ID photoionization 
models. In principle, the second of these can be tested 
against 3D photoionization code s such as MOCASSIN 
IjErcolano. Barlow fc Storevl l2005h . Until then, similar 
empirical calibrations of photoionization models across a 
range of metallicities may be the best approach. It may 
be significant that the discrepancy for the Orion Nebula 
Cluster is significantly less severe than for G29.96-0.02 
and G23.96+0.15. A further test would be for an extensive 
sample of Hn regions in the outer Milky Way and the 
Magellanic Cloud, These could potentially include 30 
Doradus (LMC) and/or NGC 346 (SMC). Beyond the 
Magellanic Clouds, spatial resolution prevents spectroscopy 
of individual stars within compact star clusters. 
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6 SUMMARY 

We present near-IR (VLT/ISAAC) and mid-IR 
(Spitzer/IRS) spectroscopy of massive stars within the 
young Milky Way cluster G10.2-0.3 (W31). Our main 
results may be summarised as follows: 

(i) H- and K-band s pectroscopy o f nake d O stars broadly 
confirms subtypes from iBlum et ail <|200lh from which a re- 
fined cluster distance (3.3 kpc) and age (~0.6 Myr) are ob- 
tained. 

(ii) W31 #26, one of the massive YSOs from lBlum et ail 
is shown to possess photospheric features, consistent 

with a subtype of 06V, in addition to near-IR circumstel- 
lar dust emission. This suggests it is in the very process of 
clearing its immediate environment. 

(iii) Mid-IR fine-structure line ratios of W31 stars overlap 
with other Milky Way and Magellanic Cloud H II regions on 
a [Ne ill] /[Neil] t o [Siv]/[Sm] di a gram. 

(iv) Following iMorisset et ail (|2004l i a comparison of 
the mid-IR radiation hardness parameter, 77(S-Ne) = 
([Nem]/[Nell])/([Slv]/[Slll] versus [Nem] /[Nell] allows de- 
pendencies upon effective temperature and ionization pa- 
rameter, U to be tested. Predicted solar metallicity (stellar 
and nebular) models differ greatly from empirical T e s and 
U, with small improvements in both for 2x solar grids. 

(v) Similar studies are planned for metal-poor environ- 
ments, which should establish whether the problem is most 
severe at high metallicity. Initial results for the ONC do 
show an improved agreement between predictions and ob- 
servations, although larger samples are required for statisti- 
cally robust results. If the discrepancy were to disappear at 
lower metallicities it would suggest that either the extreme 
UV metal line blanketing of early-types stars is incomplete, 
or there is a problem with photoionization models at high 
metallicity. 

(vi) We show that an empirical correction to the pre- 
dicted [Ne ill] /[Neil] ratio against T e g for logf/=-1.5 and 
2xZq provides a reasonable match to stellar results for W31 
and two inner Milky Way UC H II regions. Estimates of the 
O subtypes of the ionizing stars in other inner Galactic H II 
regions (Rgc < 7 kpc) are obtained. This approach is only 
practical for early- and mid- O stars, in view of the weakness 
of [Nem] 15.5/im for late O-types. 

(vii) For ground-based datasets lacking [Nem] observa- 
tions, we have used a c orrelation between [N em] /[Nell] and 
[Slv]/[Nen] - see also iGroves et ail (|2008l ) - to provide a 
calibration of [Slv]/[Nen] against T c s for super-solar com- 
pact H II regions. Since only early- and mid- O stars provide 
significant [S iv] emission, we have also obtained a similar re- 
lation for 9.0/im [Ar m]/12.8/im [Nell] versus T e s, although 
this should only be applied when other diagnostics are un- 
available, with both [Ar ill] and [Ne n] expected to be weak 
in late-type O stars. 

Finally, in view of the apparent discrepancy between 
stellar and nebular results at high metallicity compact 
Hn regions, studies of more straightforward templates are 
sought, to enable the present calibration to be put on a more 
robust footing. 
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